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Abstract 

The transient fluctuation of the prosperity of firms in a network econ- 
omy is investigated with an abstract stochastic model. The model de- 
scribes the profit which firms make when they sell materials to a firm 
which produces a product and the fixed cost expense to the firms to pro- 
duce those materials and product. The formulae for this model are parallel 
to those for population dynamics. The swinging changes in the fiuctuation 
in the transient state from the initial growth to the final steady state are 
the consequence of a topology-dependent time trial competition between 
the profitable interactions and expense. The firm in a sparse random net- 
work economy is more likely to go bankrupt than expected from the value 
of the limit of the fluctuation in the steady state, and there is a risk of 
failing to reach by far the less fluctuating steady state. 

1 Introduction 

Firms are put under the selection pressure in a network economy. Every firm 
searches for the most profitable pattern of interactions with other firms for the 
trades of goods and credits to survive. Financial distress, however, transmits 
along the interactions as well. As a firm declines in sales for whatsoever reasons, 
bad debts of the firms which sell goods to that firm tend to increase. The firms 
have a hard time to manage their financing, also decline, and some of them even 
end in bankruptcy. The bad debts increase still more. The risk that financial 
distress emerges in a network economy and transmits widely is called systemic 
risk in macro-economics. The pattern of interactions has a decisive impact on 
both the financial distress and the bankruptcy of firms. 

A supply chain and a business cluster are often found highly efficient sub- 
systems in a real network economy. A supply chain is a group of firms in a 
diverse range of industry sectors. For example, a firm in the finished goods 
sector purchases components from a firm in the component sector, and the firm 
purchases materials from a firm in the materials sector in turn. As goods are 
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moved downstream, money is transferred upstream in the opposite direction. 
The topological pattern of these interactions along a stream is a linear chain 
most simply. A business cluster is a geographic concentration of firms in closely 
related industry sectors. Empirically, firms develop reciprocally profitable inter- 
actions with each other, raise each otherfs prosperity collectively, and increase 
the productivity as a whole. Win-win games ensue there. The topological pat- 
tern of such reciprocal interactions is a ring most simply. Do chain and ring 
topologies have any special properties to keep financial distress from emerging 
and transmitting? How does the prosperity of individual firms evolve temporar- 
ily under such patterns of interactions? 

In this study, the transient fluctuation of the prosperity of flrms is investi- 
gated with an abstract stochastic model. The model describes the profit which 
firms make when they sell materials to a firm which produces a product, given a 
fixed pattern of the network for the trades of goods, and the fixed cost expense 
to the firms to produce those materials and product. The mathematical formu- 
lae of the model are parallel to those for population dynamics. The prosperity 
and the risk of bankruptcy of a firm resemble the population and the risk of 
extinction of a species. The analytic solution of the fiuctuation is presented 
for two building blocks of the entire pattern of interactions, namely, chain and 
ring with homogeneous profitability and expense parameters. The fluctuation 
and its temporal correlation are analyzed numerically for more general random 
pattern of interactions with heterogeneous profitability and expense parameters. 



2 Related works 

The contagion network model of three industry sectors describes the growth of 
firms |Delli Gatti 2010) . The model combines the change in the amount of trades 
from firms in the materials industry sector to firms in the finished goods sector, 
and the change in the interest rate of loans from banks in the finance sector to 
the firms in both industry sectors. The systemic risk emerges from the cycle that 
the cumulative loss of firms causes a sharp rise in the interest rate for the firms 
to raise funds from banks, and increases the loss still more. When the net worth 
of a firm reaches zero by netting the cumulative loss, the firm goes bankrupt. 
In the model, the bankrupt firm is replaced with a new firm. The firms and 
banks may change selectively the topology of the network for the trades of goods 
and credits as the interest rate rises and falls |Delli Gatti 2006j . The price of 
the finished goods is a stochastic exogenous variable. A large decline of price 
represents the shock which is exerted by a macro-economic environment. The 
steady state distribution of the net worth and growth rate of firms is obtained 
with simulations. These studies analyze how the avalanche of bankruptcies of 
firms ensues from exogenous macro-economic shocks, namely, a sudden fall of 
demands and a resulting decline of the price. 

Similar contagion network models are applied to analyze the transmission of 
distress among banks in an inter-bank credit network [Haldane 2011] , |Gai 2010] , 
|Martenez-Jaramillo 2010] , [Beyeker 2007 . These studies are of particular inter- 
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est worldwide since the recent financial crises. They reveal how the avalanche of 
bankruptcies of banks ensues from exogenous macro-economic shocks like a sud- 
den decrease of funds transfer requests from client firms, the collapse of an asset 
market, and the disappearing liquidity of the assets in which banks make an 
investment. None of these studies, however, addresses the problem on how the 
probability of the bankruptcy of a firm or a bank, which ignites the avalanche 
of bankruptcies, is correlated with endogenous factors like the topology of a 
network and the fluctuating trades of goods or credits between them. 

In terms of mathematical formulae, the interactions between firms are paral- 
lel to the interactions between species in the population dynamics [Allesina 2011) , 
[Larremore 20Tl] . |Wu 201 Ij . [Mehrotra 2009] . The bankruptcy of a firm when 
the prosperity reaches zero resembles the extinction of a species when the pop- 
ulation reaches zero. Jain and Krishna's linear catalytic network is a weighted 
digraph whose vertex and arc represent a species and the presence of catalysis 
from one species to another [Jain 200"T| , [Jain 1998] . The weight of an individual 
arc is the strength of the catalysis. The steady state distribution of the rela- 
tive populations of species is a function of the topology of the digraph. Species 
of the smallest population in the steady state are dismissed and replaced by 
new species. This results in an update of the topology of the digraph. An 
auto-catalytic set appears by chance after an update [Jain 2002a] , [Jain 200Tb] . 
It is a set of catalyzed species which governs the proliferation, extinction, and 
recovery of the entire system of species. Many real systems, however, do not 
reach the steady state quickly, but remain in a transient state. The strength 
of the selection pressure on individual species and the resulting probability of 
their extinction are not predicted from their relative populations in the steady 
state which these studies reveal. 

3 Stochastic model 

A stochastic model for a network economy is presented in this section. A fixed 
digraph represents the entire network economy. The vertex Vi is the i-th firm 
for i = 0,l,---,A^— 1. A'^is the number of firms. The arc Vj — > Vi represents 
the interaction from the vertex Vj to Vi. An arc Vi Vi is a loop. 

Continuous variables yi(t) is the prosperity of the i-th firm at continuous 
time t. The prosperity is an abstract quantity to represent the net worth, or 
the cumulative profit, of the firm. This is parallel to the population of a species 
in population dynamics. Assume two types of abstract interactions between 
firms. As the j-th firm produces products, its purchase of materials gives rise 
to the increase of the prosperities of a set of firms. The rate of change of j/i 
is given by Cijyj{t) where the profitability parameters Cij are constant. This 
resembles the j-th species (prey) which catalyzes the proliferation of the i-th 
species (predators). The non-zero elements of the matrix Cij represent the arcs 
E = {vj — Vi}. The decrease of iji results from the fixed cost expense by 
paying for the labor to produce products. The rate is given by —diyi{t) where 
the expense parameters di are constant. This resembles the individual deaths 
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of the i-th species in the population dynamics. 

Stochasticity ensues from an unpredictably irregular temporal pattern of 
these interactions. If the probability of an interaction per unit time is constant, 
the number of interactions in a given time interval obeys a Poisson distribution 
[Komorowski 20lTj , where the amplitude of deviation is equal to the mean. The 
time evolution of yi{t) is given by a system of Langevin equations in eq.dl]). The 
functional forms of the Gaussian white noises, ^£,j{t) and '^^i(t), are not known. 

The solutions of eq.dl]) are described equivalently by the joint probability 
density function P{y, t) of a random variable y = [yo, • • • , Vn-i) at t. Its time 
evolution is given by a Fokkcr-Planck equation [Fa 201 ij , [Dobrinevski 2012] in 
eq.jS]). 

The elements of the drift vector A and the diffusion matrix B is given by 
eq.® and ^. 

N-l N-1 

= E ^vV] = E ~ dAj)yj. (3) 
j=o j=o 

N-l 

Bij = E {(cifc + diSik)Sij + y^CikC-jkil - Sij)}yk- (4) 

fc=0 



4 Analytic solution 

A comparative analysis of the impact of topologies on the fluctuation of pros- 
perity is presented in this section. The digraph for the entire network economy 
includes a number of chain and ring sub-digraphs. A ring sub-graph is the sim- 
plest abstract form of collectively reciprocal interactions which are seen in a 
business cluster. This is an auto-catalytic set in the population dynamics. A 
chain sub-graph represents an abstract supply-chain where reciprocal interac- 
tions do not work. The nature of these building blocks is the basis to understand 
the impact of the topology of general random digraphs. The formulae for the 
fluctuation are derived for these two analytically tractable sub-digraphs with 
homogeneous profitability and expense parameters. 

The arcs of the chain are £^0 — {vq — > vi,vi — > W2, • • • ,WAr-2 ^ vn-i}- 
The vertex vq is the most downstream purchaser, and vn-i the most upstream 
seller. Assume homogeneous parameters, Cij = c > for j = i, Cij = c' > for 
j = i — l, Cij = for the other values of j, and di — d > for all i. Self-sustaining 
profitability, a firm with reciprocal subsidiaries or stably profit-making assets, is 
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Table 1: Set Q{i) which gives the bounds of summation with q and q' for ^iii{t) 
in eq-dl]). 





chain 


ring 




Q<q<i 
q' = i-q 


< g < iV- 1 
q' = i — q mod N 



meant when c > (loops {vq — > fo, • • • ,vn-i — ?> uat-i}). This is self-catalysis 
in the population dynamics. The arcs of the ring are = Eq U {wat-i — > wo}- 
Assume — c for j = i, Cij = c' for j = i — 1%N, Cij = for the others. The 
modulo operation i — 1%N outputs the remainder on dividing i — 1 hy N. 

The time evolution of the moments of y is derived from the Fokker-Planck 
equation [Maeno 2010] . [Maeno 2011| . The first order moments, ^Hi{t) = {yi)t, 
is the ensemble average of j/i at t. The initial condition is ^/ii(0) = yi{0)- It is 
given by a function of t in eq.(l5]). It increases monotonically if c > d. 

V.(0 = [exp(At)y(0)],: = J2 ^^^t^'e^-"^)*. (5) 

The bounds of summation with non-negative integer index variables q and q' 
are given by the set Q{i) in Table [TJ For the ring, q' is given by a Diophantine 
equation q' = i — q mod N whose solution is q' = i — q,i — q + N, ■■ ■ if « > g, 
and q' — i ~ q + N,i — q + 2N, ■ ■ ■ if i < q. Significant profitability from the 
nearest neighbor firm is ignited when c't ^ 1. The ignition time is tp ~ 1/c'. 
The terms of g' > iV represent reciprocal profitability. They are ignited when 
{c't)^/Nl - 1. The ignition time is t,{N) - N/ec' if TV > 1. 

The second order moments about the mean (covariance) , '^^ij{t) = {yiyj)t — 
^IJ,i{tyfj,j{t), is given by eq.®. The initial condition is ^iJ,ij{t) = 0. The operand 
functions of the three summation operators ^ are the same. 

/■* ~ ~T 
V^iW = / iexpA{t-t'){B)t'expA {t - t'))^jdt' 
Jo 

= {{c+d) J2 E +^ E } 

Qlihj) Q2(i,3) Q3(i,j)UQi(i,j) 

2/g(0)c'«'+9"+9"' 



q'lq"lq"'l 



-Tg„^,,+,,„{t). (6) 



The bounds of the ^ operators with non-negative integer index variables q, 
q', q", and q'" are given by the sets Qi{i,j) through Q4{i,j) in Tabled] Tr^r'{t) 
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Table 2: Sets Qi{i,j) through Q4{i,j) which give the bounds of the ^ operators 
with q, q', q" , and q'" for ^Hij{t) in eq.®. 





chain 


ring 




<k<i,j 


< fc < TV - 1 




0<q<k 


< g < iV- 1 


Qi{-i,j) 


q' = i — k 


q' = i — k mod N 




q" = k-q 


q" = k — q mod N 




r ^J-k 


q'" =j - k mod TV 




l<k< i,j 


< fc < TV- 1 




0<q<k-l 


< g < TV- 1 


Q2{i,j) 


q' = i — k 


q' = i — k mod TV 




q" = k-q-l 


q" = k - q - I mod TV 




r ^J-k 


q'" = i - k mod TV 




1 < k <i,j + 1 


< fc < TV - 1 




0<q<k-l 


< g < TV- 1 




q' = i — k 


q' = i — k mod TV 




q" = k-q-l 


q" = k - q - 1 mod TV 




q'" =3-k + l 


q'" = j- k + l mod TV 




< k <i,j -1 


< fc < TV - 1 




0<q<k 


< g < TV- 1 




q' = i — k 


q' = i — k mod TV 




q" ^k-q 


q" = k-q mod TV 




r - J - fc - 1 


q'" = j - fc - 1 mod TV 



6 



is specified by two non- negative integers r and r', and given by eq.(I7]). 



Jo 

^^ r'Cr"{-ir'-''"{r + r'-r")\ „ 
^ (c- 

r"=0 ^ ' 



t-\-t' — t" 



X 



r"'=0 



Eq.® and ^ hold for both the chain and ring sub-graphs. The effect of 
the reciprocal profitability of a ring appears in the bounds of summation in 
Table [1] and [2] The 3rd order moments, '^^ijk {t), and higher order moments 
are derived similarly. The explicit formula for P{y,t) is obtained by a multi- 
variate Edgeworth series [Balitskaya 1988| . The fluctuation is defined by eq.®. 
This quantity is also used in quantifying the noise in a changing environment 
[Hilfinger 2011| , and the variability of spreading phenomena [Crepey 2006| . 



Figure [1] shows Fait) for ring and chain sub-graphs as a function of t. For 
a ring, c' > alone does not contain the diverging fluctuation. Fa converges 
if c > or d = 0. Downstream vertices (vq etc.) suffer from larger fluctuation 
then upstream vertices (wat-i etc.) if d > 0. For a ring, the fluctuation in the 
transient state is larger than that in the steady state when N > 5. 

The limit for the chain is given by eq.®. It does not depend on any of i, j, 
N, c'. For any i, j, the dependence on the initial condition is ?/o(0)^^^^. Without 
self-sustaining profitability, the limit diverges unless the expense for production 
is absent. Convergence is slow t^^), and particularly slow (~ i~^) when 
c = d = 0. 

{^/^^ c-d>0 
Y yo(0) c-d 
oo c-d<0, c^O, d^O ■ (9) 

c=d=0 

The limit for the ring is given by eq. ([T0| . It does not depend on i, j. For 
any i, j, the dependence on the initial condition is ^(O))"^^^- Because 
the limit in ea.([TU| tends to decrease as N increases {N~^^'^), the gap between 
the limits for the chain and ring is more conspicuous as N increases. As c, c' 
decreases or d increases, F increases. Even without self-sustaining profitability, 
the limit is finite as far as d < c' is satisfied. 



1 c+c' +d+2Vcc' ^ _ ^ \ 

lim F,, = <( V^^^^ c a> ^^^^ 

c-d<~c' 
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Figure 1: Fii(t) as a function of t. (a) chain of = 30 (i=0 to 29) when c = 2 
(self-sustaining profitabiUty), c' — 1, and d — 1. For all i, yi{Q) = 10. (b) chain 
when c = 0, c' = 2, and d = (no expense). Nothing happens for i ~ Q. (c) 
chain when c = 0, c' = 3, and d = 1. (d) rings of TV = 2, 3, 5, 10, 20, 30, 50, 
and 100 when c = 0, c' = 3, and d ~ 1.^ 



At small t, Fij is given by ea. ([TT|) . The time evolution of Fij for chains is 
equal to that of rings of any length with the same parameters. Vertices are 
exposed to nearly identical trends of growing fluctuation. None of these sub- 
digraphs are discriminated from each other in the digraph with homogeneous 
parameters in the early transient state. 



5 Numerical analysis 

The nature of more general random digraphs is presented in this section. First, 
the fluctuation for digraphs with heterogeneous profltability and expense pa- 
rameters is analyzed by evaluating the integral in eq.® numerically. Fig. [2] 
shows Fii[t) for a chain, a ring, and a chain of rings when is chosen ran- 
domly between (1 — h)c and (1 -f ft,)c, given the median c and heterogeneity 
parameter h. Given d and y(0), di and j/i(0) are chosen randomly in the same 
manner. 

In (a), most curves for F take a local maximum at about t — tp which is fol- 
lowed by a few local extrema. The number of the extrema depends on vertices. 
The curve for an downstream vertex (with i) along the chain tends to attract 
those for its upstream neighbor vertex (with i + 1 etc.). Fu < F^k holds for any 
l> k bX large t. The degenerate limit in eq.® is split into multiple values. The 
number of the limits depends on the topology and the value of the parameters. 
Similar numerical examples show if the chains without self-sustaining profltabil- 
ity are shorter, the expense starts surpassing the profitability and an infinitely 
growing trend of F appears earlier. In (b), the curves for F for the ring oscillate 
regularly after about t = t,-. Crests and troughs of the profltability arise locally, 
and transmit along the ring as a travelling wave. This profit-making wave does 
not impair the convergence to the single limit. Generally, if the rings are shorter, 
the reciprocal profltability is ignited and F starts converging earlier. In (c), the 
fluctuation resembles that for the chain, but the number of limits tends to be 
smaller. The swinging changes in F at about t — tp until about t represent 
the transient state from the common initial growth to the flnal steady state. 

Next, the dependence of the transient state on the number of arcs and 
rewiring of random digraphs, and on the heterogeneity of profltability and ex- 
pense parameters is investigated with the temporal correlation of Fait). The 
correlation is the Pearson's product- moment correlation coefficient of F between 
the inception of growth at ti = 0.01 and t. It is defined by ea. ([T^ . F{t) is the 
average of Fu over the all vertices at t. 




+ 



(11) 



j:^iFu{t)^FmFu{h)-F{h)) 



(12) 



j:AFu{t)~F{twj:^{Fu{h)~F{h)r 
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(a) 0.6 




' ' — ' — ' — • — 

0.01 Q.I 1 10 100 




0.01 0.1 1 10 100 




Figure 2: Fii{t) for z = to 29 {N = 30) as a function of t when aj, di, and 
j/i(0) are chosen randomly, (a) Chain with c = 2 for j = i (self-sustaining 
profitability), 1 for j = i — 1, and for the other j. (b) Ring with c = 3 for 
j = i — 1%N and for the other j. (c) Chain of rings whose topology is shown 
in (d) with c = 3 if an arc is present and otherwise. For the all graphs, d = I, 
y{0) = 10, and h = 0.5. 
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Figure 3: R{t) as a function of t for random digraphs of = 100. (a.) K — 2 
and h = 0.05 with c = 3 if an arc is present and otherwise, (b) K = 2 and 
h = 0.05 with c = 2 for j = I (self-sustaining profitabihty), 1 if an arc is present 
and otherwise, (c) Ring oi K = 5 and 20 with c = 3 if an arc is present and 
otherwise. For the all graphs, d = 1 and y{0) = 10. The rings with = 1 are 
Erdos-Renyi random digraphs. 
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Fig. [3] shows R{t). The curves are the average over many different digraphs 
with the same number of arcs, rewiring, and heterogeneity of parameters. The 
graph gives the general tendency of the fluctuation as a function of K and 
given h. A digraph is specified by K and Pr- K is the out-degree of vertices 
which is the average over the all vertices. Sparse digraphs have small values of 
K. Some of the arcs are chosen at the rewiring probability of p^- Either end 
of the chosen arcs is rewired to a different vertex. The figure shows the curves 
for rings with pr = 0, 0.01, 0.1, and 1, namely, random digraphs in the Watts- 
Strogatz model [Watts 1998] . and chains with p^ ~ 0. The Watts-Strogatz 
model describes the small-world property of digraphs. The rings with = 1 
are equivalent to random digraphs in the Erdos-Renyi model |Erd6s 1959] . In 
the Erdos-Renyi model, an arc is set between each pair of vertices with equal 
probability, independently of the other arcs. 

In (a), the sign of R turns about sharply after the ignition of profitability. For 
the rings with small Pr, the initial correlation (R ^ 1) is followed by moderate 
anti-correlation {R > —0.5), which disappears as the reciprocal profitability 
takes over the role to govern F from the profitability and expense. The anti- 
correlation decreases only very slowly for the rings with Pr = 1 where rewiring 
destroys collectively reciprocal interactions. Strong anti-correlation {R ^ —1) 
appears for the chains where the profitability and expense keep on working for 
ever. Both the profitability and expense cause anti-correlation, and the anti- 
correlation is a potential indicator of the transient state. The reason for these 
is as follows. Initially, a small prosperity means relatively large fluctuation 
because of the factor ^y {dyi{0) + c'yi_i(0))/j/i(0) in cq.lfTTj) when c = 0. After 
the ignition of the profitability, it mitigates large F for the vertex having small 
y(0). The expense pressures small F for the vertex having large y{0) into rising. 
So the anti-correlation appears. At the same time, they tend to compete with 
each other in forcing the absolute value of the individual Fa to decrease or 
increase. 

In (b), anti-correlation disappears gradually as a result of self-sustaining 
profitability for the chains. This effect is not evident for the rings. In (c), 
if the values of K and h decrease, or if the value of Pr increases, the sign of 
R turns about earlier and subsequent anti-correlation becomes stronger. The 
time when the correlation converges also depends on them. The fluctuation 
for sparse random digraphs undergoes an upheaval in the salient and long-lived 
transient state because of the absence of collectively reciprocal interactions. 
When Pr = 1, the probability at which at least one path starting from a vertex 
forms a ring sub-digraph of length L (< log^ iV) is given roughly by Pa.{L) k, 
1 - {{N -2)/{N - 1)}^"" if TV > 1. When N = 100 and K ^2, p^{2) w 0.04. 
Ring sub-digraphs are missing in the digraph. When K = b, Pa,{2) « 0.22. Ring 
sub-digraphs appear. When K = 20, Pa{2) « 0.98. Most vertices in the digraph 
are vertices of ring sub-digraphs. If K is larger, shorter ring sub-digraphs appear 
more and the convergence of F starts earlier. 
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6 Discussion 



The fluctuation does not converge to the Umit monotonically. If a firm is a 

more downstream purchaser (with smaUcr ?') along a siipply-chain, the expense 
starts surpassing tlie profitabihty tliere earlier. The supply-chain would not be 
stable for ever without self-sustaining profitability. This means that there is a 
finite life time for supply-chains. If the business cluster to which a firm belongs 
is smaller, the reciprocal profitability is ignited earlier. Just the presence of a 
business cluster does not predict correctly the risk of impending bankruptcy of 
inclivichial firms in the transient state. A mimber of such biiilcling blocks are 
inter-connected with each other to form a general digraph for the entire net- 
work economy of firms. The time when the fluctuation converges depends on 
the sparscness of the general digraph. The swinging changes in F arc the conse- 
quence of a topology-dependent time trial competition between the profitability 
and expense. 

No doubt it is a rational strategy to survive in the long run for any firms to 
occupy a position in a predominant business cluster with reciprocal interactions, 
say, an intersection of multiple long ring sub-digraphs in the digraph of the entire 
network economy. But in practice, every time the firms settle accounts, their 
prosperity must be positive. The selection pressure is brought on the flrms 
nearly continuously if the time interval between their successive settlements 
is short. Even if a business cluster is present topologically, it may not take 
effects, nor will it start working soon. The flrms must comprehend the time 
trial competition between the profitability and expense, and determine whether 
they remain at the current position or bet on the success of rewiring to push 
back against the immediate selection pressure. 

The bankruptcy of a firm ensues directly from an instantaneous fluctuation. 
It is not within the scope of this study to solve the problem of a flrst passage time 
with reflective boundary conditions at yi = 0, but qualitatively, the bankruptcy 
of the i-th flrm is more probable as Fa increases. If P(y,t) is approximately 
a multi-variate normal distribution, the rough estimate of the probability of 
bankruptcy at t is ph{t) ^ i_^PM{yi,t)<lyi = ^{-l/Fu{t)) where Pm is a 
marginal distribution of P and $ is the cumulative density function of a normal 
distribution. F in the transient state often exceeds the limit of F. The small 
value of the limit of F^ does not mean the certainty of the survival of the i-th 
firm until then. Some firms are more likely to go bankrupt than expected from 
the fluctuation in the steady state, and consequently, there is a risk of failing 
to reach by far the less fluctuating steady state. The risk is more significant in 
sparse random digraphs. 
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